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Abstract. We detail the experimental observation of the non-equilibrium many-body 
phenomenon prethermalization. We study the dynamics of a rapidly and coherently 
split one-dimensional Bose gas. An analysis based on the use of full quantum 
mechanical probability distributions of matter wave interference contrast reveals that 
the system evolves towards a quasi-steady state. This state, which can be characterized 
by an effective temperature, is not the final thermal equilibrium state. We compare 
the evolution of the system to an integrable Tomonaga-Luttinger liquid model, and 
show that the system dephases to a prethermalized state rather than undergoing 
thermalization towards a final thermal equilibrium state. 
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1. Introduction 

A general understanding of the dynamics of non-equilibrium many-body quantum 
systems is an important unsolved problem impacting on many areas of physics [1, 2]. 
In particular the question of why and how isolated quantum systems relax towards 
equilibrium states has only been studied for a very limited number of special systems [3, 
4, 5] . Moreover, in many cases, such as when dealing with integrable systems exhibiting 
many constants of motion, efficient relaxation can be completely absent [6] or strongly 
inhibited. 

An intriguing phenomenon which has recently been experimentally demonstrated 
in this context is prethermalization [7, 8]. The concept of prethermalization was first 
introduced in order to explain the remarkable success in applying thermodynamic 
models to heavy-ion collision experiments on time scales much shorter than the 
expected thermalization time of the scattered particles [7, 9, 10, 11]. In particular, 
the successful application of hydrodynamics [12] relies on the assumption of being at 
least approximately close to thermal equilibrium at every point in space and the presence 
of an equation of state [13, 14]. 

Taking a low-energy quark-meson model as a model system, Berges et al. [7] 
suggested in 2004 an explanation for this unexpected, apparent early thermalization. 
In their theoretical analysis, Berges et al. were able to show that after a very rapid 
initial evolution, their model system relaxes to a non-thermal quasi-steady state. This 
intermediate state is very robust with respect to changing initial conditions and already 
shows many of the bulk characteristics of the final thermal state. For example, the 
ratio of kinetic energy to pressure (i.e. the equation of state) or a properly defined 
kinetic temperature are already practically indistinguishable from their values in thermal 
equilibrium. In contrast, mode quantities such as the occupation numbers of different 
momentum modes are still far away from a true thermal equilibrium even though during 
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the initial rapid evolution they also reach quasi-steady values. Their work showed that 
although the model system had not yet fully thermalized, it could in many aspects be 
described using a thermal model. Accordingly, the above described phenomenon was 
termed pre-thermalization. 

Since then, prethermalization has been discussed in the context of many other 
systems, ranging from models in relativistic field theories [15], the dynamics of the 
early universe [16] to systems which can be realized using ultracold atoms, such as 
systems obeying the Hubbard model [17, 18, 19], spinor condensates [20], quantum 
Ising chains [21], two-dimensional superfiuids [22] or one-dimensional (Id) quasi- 
condensates [23, 24]. 

In the present understanding, prethermalization is characterized by the rapid 
establishment of a quasi-stationary state on time scales much shorter than the expected 
thermalization time. This state is long-lived and already exhibits some thermal-like 
properties but can still be very different from the true thermal equilibrium state of the 
system. Full relaxation to thermal equilibrium, if present at all, is then expected to 
happen on a further much-longer time scale. While thermalization can be interpreted 
as a full loss of the memory of the system about its initial state, prethermalization 
describes only a partial loss of information about the initial state. 

In the following, we discuss the direct observation of the prethermalized state 
presented in [8], giving, particularly, details of the experimental system and analysis 
methods used to reveal the prethermalization [25]. 

2. Non-Equilibrium Dynamics of Interfering One-Dimensional Bose Gases 

2.1. Experimental System 

The experimental study is performed using trapped Id Bose gases. Such systems offer 
two unique advantages for non-equilibrium experiments: firstly, on the experimental 
side, realizing them with ultracold atoms facilitates a precise preparation and probing 
of the system. Secondly, on the theoretical side. Id Bose gases offer a model system 
which contains complex many-body physics but can still be captured with reasonable 
theoretical effort, particularly due to the existence of effective models which allow 
to describe the essential physics in a relatively simple way [26]. Furthermore, the 
homogeneous Id Bose gas with repulsive contact interactions is an example of a fully 
integrable quantum system [27, 28]. The approximate realization of such a system 
in experiments thus allows the study of thermalization in the vicinity of multiple 
conserved quantities and hence the study of the interplay between integrability, many- 
body dynamics and thermalization. 

Figure 1 summarizes the main idea of our experimental study of prethermalization 
using Id Bose gases in the quasi-condensate regime [29, 30, 31]. In this regime, 
density fluctuations are strongly suppressed and the gas is characterized by strong 
phase fluctuations. The properties of these phase fluctuations are determined by the 
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temperature and the density of the system. In the experiment such a single trapped 
quasi-condensate is rapidly and coherently split, producing two uncoupled Id gases 
with identical phase profiles. The aim of our study is to probe how these initially 
almost perfect correlations of the relative phase become obscured over time and if the 
thermal equilibrium state corresponding to two completely independent gases is finally 
reached [32, 33]. To this end, the two gases are allowed to evolve in the double- 
well potential for a varying evolution time before the relative phase correlations 
are probed via time-of- flight matter wave interference (figure Ic). As differences in 
the relative phase lead to a locally displaced interference pattern, the contrast of the 
longitudinally integrated interference pattern is a direct probe for the dynamics of the 
system [34, 35, 36, 37]. Example interference patterns after various evolution times, 
demonstrating the loss of the initial phase coherence, are shown in figure 1. 
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Figure 1. Overview of the experiments performed in the context of this work [8]. 
(a,b) A phase fluctuating Id quasi-condensate is coherently split creating two Id gases 
with almost identical phase profiles 4>l{z) and 4>r{z) (represented by the black solid 
lines). The gases are then allowed to evolve in the double- well potential for some time 
which leads to strong fluctuations in the local phase difference 0(z) and a decrease of 
the phase correlation length A^. The question at which this experiment aims is whether 
and how this dynamical state reaches the thermal equilibrium state corresponding to 
two independently created quasi-condensates. In this thermal equilibrium state the 
phase difference between the Id gases fluctuates strongly along its length and the 
correlation length At is determined by the temperature and density of each cloud, (c) 
The phase difference 4){z) between the two Id gases is probed through time-of-flight 
matter wave interference of the two gases where the local relative phase is directly 
transformed to a local phase shift of the interference pattern. The contrast C{L) of 
the axially integrated interference pattern can then be used as a direct measure of the 
strength of the relative phase fluctuations. 
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2.2. Initial N on- Equilibrium State and Thermal Equilibrium 

A prerequisite for non-equilibrium experiments is the ability to precisely prepare and 
characterize both the initial non-equilibrium state, as well as the expected thermal 
equilibrium state of the system. One of the key advantages of coherently split Id Bose 
gases is that both these states can be prepared and described with high precision. 

For a general system of two spatially separated Id Bose gases, the excitations can 
be described by anti-symmetric and symmetric longitudinal modes which relate to the 
relative and common degrees of freedom of the two halves of the system. They are given 
by 

m = My)-My) , <l>.oM = "^^^^^^^ (i) 

for the phase, and 

n{y) = ^^^^^ ^ ^^^^^ , Ucomiy) = riniy) +nL{y) (2) 

for the density. Here 4>L,R{y) describes the longitudinal phase profiles and rii jiiy) the 
density of the left and the right gas, respectively. 

After the coherent splitting, the individual phase profiles of the two halves are 
almost identical and hence the relative phase profile is almost fiat (figure 2 left column). 
In terms of excitations, this means that all the thermal excitations are initially contained 
in the common degrees of freedom. As we will detail in section 6, the relative degrees 
of freedom, on the other hand, are initially populated only by quantum noise created in 
the splitting process. 

In thermal equilibrium we expect to find the system in a state where the energies 
in the relative and in the common degrees of freedom are equal (figure 2 right column). 
This thermal equilibrium state corresponds to the situation of two independently created 
quasi-condensates where the phase profiles of the clouds are uncorrelated down to the 
thermal phase correlation lengths, resulting in a relative phase profile that fluctuates 
strongly along the length of the system [32, 33]. In experiment, the thermal equilibrium 
situation can be purposely created, by splitting a thermal gas into two, followed by the 
creation of two independent quasi-condensates through further cooling in the double- 
well. 

The transition between the initial non-equilibrium state and the thermal 
equilibrium state can therefore be directly probed by studying the relative phase 
fluctuations. 

Considering above, we note that any thermalization mechanism must necessarily 
redistribute the mode populations of the relative and common degrees of freedom to 
lead to the thermal equilibrium situation. 

2. 3. Full Distribution Functions of Matter Wave Interference Contrast 

Previous experiments have also investigated the dynamics of split Bose gases [38, 39] 
or equivalent systems [40]. In these works, though, due to the lack of theoretical 
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Figure 2. Comparison of the initial non-equilibrium state prepared by the coherent 
splitting process and the thermal equilibrium state of two independent quasi- 
condensates. The system can be described by a combination of relative degrees of 
freedom (a) and common degrees of freedom (b). After the coherent splitting the 
common degrees of freedom contain all the thermal energy E ~ fc^T of the system, 
whereas the relative degrees of freedom contain only the quantum noise of the splitting 
process. Here, kg is Boltzmann's constant and T the temperature of the system. In 
thermal equilibrium all degrees of freedom contain the same amount of thermal energy 
and are thus characterized by the same temperature. 



or experimental tools for a full characterization of the transient states, the driving 
mechanism of the dynamics could not or only indirectly be revealed. Consequently, the 
state to which the system decayed remained elusive. 

This directly illustrates the key challenge common for most experiments trying to 
observe non-equilibrium dynamics: the scarcity of experimental tools for characterizing 
the complex many-body states during the evolution. In this work, this problem is 
approached by generalizing tools developed for studying equilibrium systems of spatially 
separated Id Bose gases [35, 32, 41, 33] by applying them to the non-equilibrium case, 
as recently proposed in [23, 24]. For this purpose, we measure the time evolution of full 
quantum mechanical probability distribution functions (FDFs). 

When dealing with many-body quantum systems, valuable information about 
the underlying quantum state of the system can be gained through noise correlation 
measurements. Performing such measurements already deepened our understanding of 
quantum mechanics as it led to the discovery of the Hanbury Brown-Twiss effect [42] 
which triggered the development of modern quantum optics [43]. Furthermore, the study 
of current fluctuations led to important observations in quantum-Hall systems [44, 45] . 
Recently in atomic physics the analysis of noise correlations revealed the coherence 
properties of atom lasers [46] and enabled observations of the Hanbury Brown-Twiss 
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effect for massive fermions and bosons [47] . It was furtfier suggested [48] and 
experimentally demonstrated [49, 50, 51, 52, 53] that noise correlations in time-of-fiight 
can be used to probe strongly-correlated equilibrium states of many-body quantum 
systems. 

In Id systems fluctuations play a much more pronounced role than in their three 
dimensional (3d) counterparts [54, 55]. The interference pattern of two expanding 
Id quasi-condensates therefore inherently contains strong noise and fluctuations that 
can be used to study the many-body state of the system. When performing 
interference experiments with spatially separated Id Bose gases, one is furthermore 
in the advantageous position that the interference pattern itself is directly related 
to the correlation functions of the system. For example, the mean squared value of 
the integrated matter wave interference contrast (|Cp) is a measure of the integrated 
two-point phase correlation function. Similarly, by exploiting the strong shot to shot 
fluctuations of the contrast, higher moments of the contrast can be obtained which 
contain information about higher-order correlation functions [35, 56]. An experimentally 
accessible quantity capturing all these higher moments of the contrast and therefore also 
the higher-order correlation functions is given by the full distribution function of the 
contrast P{C). This function describes all the fluctuations of C in the measurement, 
as -P(C) dC measures the probability of observing a contrast between C and C + dC. 
Higher moments of {C^) of the contrast can be obtained by integration j C^P{C) dC. 

It is important to point out that the method of using FDFs to characterize a system 
requires the detection of single realizations of the quantum system in question. If only 
ensemble averages can be measured in the experiment, the statistics of those values 
will always be Gaussian due to the central limit theorem and the characteristic higher 
moments of the observable will not be accessible. The method of analyzing quantum 
states through the FDFs has already been successfully used to characterize isolated [32] 
and tunnel-coupled [33] systems of two separated Id Bose gases in equilibrium. Within 
the context of [23, 24], this method was extended theoretically to non-equilibrium 
systems, whereas the experimental demonstration of this method is the topic of the 
work presented here and in [8]. 

3. Experimental Details 

The preparation of the atomic sample is performed following our standard procedure [57] 
using magnetically trapped ^^Rb in the 5S'i/2 F = 2, nip = 2 state on an atom 
chip [58]. The atom chip is a current-carrying gold structure micro-fabricated on a 
silicon substrate. A basic outline of the most important structures on the chip is given 
in figure 3a. 
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Figure 3. (a) Schematic of the atom chip layout. (A) The 100 /im-wide central 
trapping wire, (B,C) 30 /im-wide wires for RF currents, (D) the atom cloud, (E) the 
silicon substrate. Additional structures, at a distance of ~ 1 mm from the condensate, 
are used to provide longitudinal confinement, (b) Simulation of the longitudinal density 
profiles of trapped quasi-condensates gases following [59]. The black solid line shows 
the profile of the cloud before splitting for an atom number of 6000 atoms. The colored 
dotted, dotted dashed and dashed lines show the profile of the sum density of two quasi- 
condensates trapped in the double-well trap with a longitudinal trapping frequency of 
6.5, 7, and 7.5 Hz. For the experimentally realized double-well trapping frequency of 
7 Hz, the profile matches almost exactly to that of the single- well trap at 11 Hz. 



3.1. Realizing a Quench: The Coherent Splitting Process 

We first create a single-well trap using the static magnetic field of a 100 /xm wide wire 
in combination with homogeneous bias fields. Longitudinal confinement is provided by 
additional wires. The trap is located at a distance of about 100 /im away from the atom 
chip surface (see figure 3a). This strongly anisotropic trap has measured trap frequencies 
of u;_L = 27r X (2.1 ± 0.1) kHz (radial confinement) and u\i^ = 27r x (11 ± 1) Hz (axial 
confinement). 

To create the double-well potential we use an adiabatic dressed-state potential [60, 
61]. By applying RF current to two 30 yum-wide wires that are adjacent to the central 
100 /xm wire, the cloud is split along the radial direction of the trap and perpendicular 
to the direction of gravity. This ensures that any potential sag due to gravity is common 
for both wells and that the 1/r dependence of the RF fields emanating from the wires 
is the same for both wells. For the splitting process, the amplitude of the RF current 
is linearly increased from to typically 22.5 mA within 17ms and the frequency is 
30 kHz detuned to the red of the F = 2,mF = 2 ^ F = 2,mF = 1 transition at the 
minimum of the initial static magnetic trap. This creates a double- well with a separation 
of (2.75 ± 0.05) /im, a simulated barrier height of (2.9 ± 0.1) kHz and measured trap 
frequencies for each well of u± = 27t x (1.4 ± 0.1) kHz and u\\ = 2it x (7 ± 1) Hz. Great 
care was taken to guarantee stability and repeatability of the sphtting process, especially 
to gain a symmetric splitting with almost zero mean atom number difference between 
the two halves. For the fine tuning of this quantity the relative amplitude of the RF 
current applied to the two 30 /xm-wide wires can be slightly adjusted. 
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Following [33, 62], any residual tunnel coupling between the two halves of the system 
was calculated to be J <^ 27r x 0.1 Hz. In this regime, a residual tunnel coupling would 
have no significant effect on the results presented in this work [41]. The absence of tunnel 
coupling was also confirmed by experiments with independently created condensates in 
the same double- well [33] . In the measurements presented in section 5 we further used a 
larger final RF amplitude which resulted in a larger splitting distance and larger barrier 
height corresponding to an even lower bound for the residual coupling and found no 
difference to the previous results. 

Although the splitting RF current was ramped up in a time of typically 17 ms, 
the actual splitting process is much faster than the timescale of the longitudinal 
dynamics and happens close to the end of the ramp when the coupling between the 
two gases vanishes. From the evolution of the mean relative phase between the two 
gases for deliberately imbalanced double- wells [63], we determined this decoupling to 
occur (15 ± 0.5) ms after the start of the RF ramp. All evolution times given in the 
following refer to this point in time. To capture also the evolution that happens after 
the sphtting but still during the ramp up, it is possible to probe the atoms during the 
ramp up. 

While the actual splitting is much faster than the ramp up of the RF current, the 
trap is nevertheless continuously deformed throughout the whole ramp up time. The 
choice of 17 ms total duration is a compromise between realizing a fast splitting and 
avoiding the excitation of collective oscillations of the cloud. The trapping frequencies 
of the final double-well were chosen to be such that the longitudinal profile of the 
clouds before and after splitting are almost exactly matched, as can be seen from the 
simulations in figure 3. This minimizes the creation of longitudinal breathing oscillations 
of the clouds during splitting. As we will discuss in section 7 these breathing oscillations 
and the general instability of the splitting process severely limited previous studies of 
coherently split Id Bose gases [38]. We find that longitudinal center-of-mass oscillations 
of the clouds after sphtting can not be completely avoided without the use of more 
complex splitting protocols. However, these small oscillations do not significantly alter 
the physics as they do not change any significant parameters of the cloud such as density 
or total size. 

3.2. Probing the System 

The system is studied using absorption imaging, which integrates the 3d structure of 
the expanding atomic density distribution along one direction of space. Figure 4 shows 
schematic views of the integrated density distribution as seen from each of the three 
detection directions realized in the experiment. 

Imaging along the x-direction is used to probe the system from its transversal 
direction. This can be applied both to extract information about a single quasi- 
condensate before the splitting process as well as about the incoherent sum of two 
quasi-condensates expanding from the double-well. As the phase of a quasi-condensate 
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fluctuates in space, density fluctuations will appear along its length when it is released 
from the trap and allowed to expand in time-of-flight [64, 65, 66]. These spatial density 
fluctuations are a direct consequence of matter wave interference between the different 
positions along the single quasi-condensate which have different phases. In equilibrium, 
the in-situ phase profiles are determined by the temperature of the system, so that the 
measurement of the spectrum of density fluctuations in time of flight can be used for 
thermometry of the system before the splitting process, as experimentally demonstrated 
in [66]. All initial temperatures mentioned in the following, as well as the temperatures 
of pairs of independently created gases in the double-well have been determined using 
this method. 

Information about the relative phase between the two Id clouds can be obtained 
from absorption images along the longitudinal axis of the two gases in the y-direction, 
or transversally along the vertical z-axis. 

In an intuitive picture, the interference pattern can be thought of as an array of 
many thin interference patterns stacked up together along the longitudinal direction 
of the cloud, where the phase of each thin interference pattern depends on the local 
in situ phase difference between the clouds. The axially integrated contrast C{L) of 
the interference is therefore a direct measure for the relative phase fluctuations within 
the integration length L [34, 35, 36, 37]. In the following we will discuss the technical 
implementation of the detection systems for the longitudinal and vertical directions in 
detail and show how we extract the integrated contrast C{L) for different integration 
lengths L. 

3.2.1. Optical Slicing Method In this scheme, the imaging beam propagates along the 
longitudinal axis of the Id gases and directly performs an integration of the interference 
pattern along the length of the gases. In order to vary the integration length L we 
perform the following slicing scheme using optical pumping (figure 5a). Shortly (1ms) 
after the two halves of the split condensate are released from the trap, a shadow of 
width L is imaged on the atoms using a 50 — 125 /is pulse of optical pumping light. 
The frequency and the polarization of the light are chosen such that the atoms, initially 
being in the 55*1/2 -F = 2, = 2 state, are excited to a state in the 5P3/2-F = 1 manifold 
(figure 5b). From there, the atoms preferentially decay to a state in the 55*1/2 = 1 
manifold which is dark to the imaging beam operating on the 55*1/2 -F = 2 — 5P3/2-F = 3 
transition. Atoms which are in the protected shadow area of length L are not optically 
pumped and can be imaged after whichever chosen time-of-flight using an absorption 
imaging system operating along the axial direction of the system. 

The branching ratios of the 55i/2 -F = 2 — )■ 5P3/2 F = 1 pumping transitions are 
such that, on average, only very few photons are needed for pumping. Already after 
one photon has been absorbed and remitted an atom has a probability above 80% of 
being in a dark state. This ensures that the atoms remain virtually undisturbed by the 
optical pumping process. 

To create the sharp-edged shadow for the optical slice, we use a multi-lens high- 
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Figure 4. Schematic view of the experimental situation considered. Sphtting of the 
gases is performed along the x-axis, the y-axis corresponds to the longitudinal direction 
of the Id Bose gases, gravity points downwards along the z-axis. The atomic density 
distribution in time-of-flight can be imaged along each axis of the coordinate system 
using absorption imaging. Imaging along the a;-axis is used to transversally probe single 
gases before the splitting or the incoherent sum of the two independent gases in the 
double well, the latter case being depicted in the figure. The corresponding absorption 
picture shows strong density fluctuations which are used for thermometry (details 
see text). Imaging along the y-axis and the z-axis is used to probe the interference 
pattern of two Id Bose gases. Imaging along the y-axis results in a projection of the 
interference pattern along its length, while imaging along the z-axis reveals the full 
undulating structure of the interference pattern. 



resolution imaging objective in reverse. The imaging system has a maximum resolution 
of 2.5 /xm, and a demagnification of approximately 1:4.6. The shadow is produced by a 
glass target with a series of opaque bars (produced at the ZMNS TU-Wien). In order to 
be able to change the integration length L in each experiment the bar target is mounted 
on a computer-controlled motorized translation stage. 

The effect of the optical pumping light on the detected atomic density is presented 
in figures 5c and 5d. Figure 5c shows examples of sliced atomic density distributions 
after optical pumping and a subsequent time-of-flight. The images were obtained using 
the same imaging system (x-direction) that is used in reverse for the optical pumping. 
Example interference images are shown in figure 5d for three integration lengths and a 
hold time of 7 ms in the double- well potential. Note that due to the evolution in the 
double-well potential, the contrast of the interference pattern which is integrated over 
the full cloud is already significantly reduced in most of the images. Furthermore, since 
for shorter L more and more atoms are cut away, the detected atom number becomes 
lower, giving a limit to the shortest integration length which can be investigated. 
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Figure 5. (a) Illustration of the working principle of the optical slicing method. The 
atom clouds are released from their trap and allowed to fall under gravity. Shortly 
after release, the sheer beam containing a shadow of length L is projected onto the 
atom clouds to pump the unwanted edges into a state that is dark to the imaging beam 
light. Standard absorption imaging is then used at a much later time of fight to image 
the remaining undisturbed part of the cloud, (b) The relevant transitions of the D2 
line of ^^Rb for optical pumping. Numbers indicate relative transition strengths using 
the normalization from [67] . (c) Images of the effect of slicing on a quasi-condensate as 
seen from the direction of the slicing beam, (d) Interference images and corresponding 
line profiles for different integration lengths L after an evolution time of 7 ms in the 
double-well potential. 

3.2.2. Direct Imaging Method In the vertical direction we can directly image the 
undulating structure of the interference pattern and freely select the integration length 
via post-processing. To this end an imaging system is employed which observes the 
atoms from below (figure 6). 

The implementation of such an imaging system is challenging since the chip surface 
blocks the line of sight, hindering standard absorption imaging schemes. One possibility 
to circumvent this problem is to use a light sheet fluorescence detector as in [33, 68]. 
Here, we instead rely on absorption imaging, as it allows for a higher imaging resolution. 

To this end, we reflect the imaging beam from the gold surface of the atom chip 
close to normal incidence. In this procedure the high quality of the gold layer is essential 
to achieve good images. For a collimated imaging beam the reflection means that the 
light passes the atom cloud twice [69]. This results in a second, virtual image behind the 
atom chip, which is out of focus and overlaps with the primary image. Since the imaging 
of the interference pattern requires the imaging beam to propagate almost parallel to 
the fringes, it is not possible to offset the two images by using a large angle of incidence 
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Figure 6. (a) Schematic of tlie direct imaging method, (b) Exemplary image. As 
the incident beam is focused, only a very small area around the atoms is illuminated 
(like a spotlight), whereas the rest of the CCD image is filled with unimportant noise. 
This removes a second, out of focus image which would be created if the incident beam 
passed the atoms twice. 



when propagating the imaging beam onto the atom chip surface. In order to ehminate 
the second image we therefore focus the imaging beam close to the atoms, avoiding the 
first absorption process as shown in figure 6a. This hmits the field of view to the area 
of a 'spotlight' (figure 6b) and ensures that the light only interacts with the atoms once. 
Also, by aligning the slowly diverging beam onto the central 100 /im wire of the chip, 
interference effects caused by diffraction from the microscopic chip structures can be 
minimized. We note that a similar method was used in [70]. In figure 6b we show an 
experimentally obtained image, where the imaging system was set up for a time-of-flight 
of 16 ms. The image clearly reveals the full structure of the spatially varying interference 
pattern, which can then be used for further analysis. 



3. 3. Extraction of the Matter Wave Interference Contrast 

The main observable used in this work is the contrast C{L) of the matter wave 
interference pattern integrated over a length L (figure Ic). For the optical slicing the 
integration length L is selected in the imaging setup by choosing a specific mask, and 
the integration is performed optically. In the direct imaging method the recorded image 
of the interference pattern is integrated over a the length L during post-processing. In 
both cases this results in a line profile (see figure 1). This line profile is fitted with a 
sinusoidal function fhix) having a Gaussian envelope: 



fdx) 



1 + c{L) cos ( ^ + e, 



(3) 



A 

where u is the cloud size, A is the fringe spacing and 6'^, is the global phase of the 
interference pattern integrated over L, i.e. the phase of the line profile. As interactions 
do not play a significant role in the radial time-of-flight expansion for our parameters, 
the fringe spacing A can be directly related to the double- well-potential separation 
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d = ht/rriA, where t is the time-of-flight, m is the mass of the atoms and h is Planck's 
constant [71]. Thus, the stabihty of the double- well splitting can be monitored with 
interferometric precision. The contrast C{L) is extracted from this fit and used for 
further analysis that is presented in the remainder of this paper. Repeated realizations 
of the experimental cycle then allow us to build the time-dependent FDFs of the 
interference contrast C{L). 

3.4- Comparison of the Two Imaging Methods 

The use of two technically different imaging systems allows us to preclude systematic 
effects of the imaging system on the observed contrast evolution. 

Example FDFs obtained for different integration lengths L from both imaging 
systems are shown in figure 7. The total area under the curves is normalized to 1 
in order to represent a probability density. The upper panel shows FDFs of the squared 
contrast C^. While the overall form of the FDFs show qualitative agreement, a small 
offset between the results of the two imaging systems remains. This is due to the fact 
that the absolute value of is affected by the finite imaging resolution and any mis- 
alignment of the imaging system. These effects can be eliminated by analyzing contrast 
distributions that are normalized to the mean contrast (C^). After this procedure the 
offset vanishes, as shown in the lower panel of figure 7. 

This comparison confirms that the form of the C^/ (C^) distributions is not 
influenced by any technical effects of the imaging systems and can therefore be used 
to study the dynamics of the system. Note that the normalization to the mean contrast 
also removes further systematic errors affecting both imaging systems, as for example 
imperfections in the switch-off of the trapping fields. 

In the following, the only and important practical difference between the two 
methods is thus that the total measurement time with the direct imaging is significantly 
shorter than the time needed to obtain the data with the optical slicing method. This is 
a consequence of the fact that with the direct imaging method only one experimental run 
is necessary to extract information on all length scales down to the imaging resolution. 
Systematic studies of large parameters spaces or long evolution times, as presented in 
section 5, are thus always performed with the direct imaging system. 

4. Observation of Prethermalization 

4.I. Dynamics of the FDFs of Matter Wave Interference Contrast 

Figure 8a shows the measured time dependence of the FDFs of interference contrast 
C{L) for the double- well system which was initialized through a coherent splitting 
process, as described in the previous section. The data was obtained using the optical 
slicing detection scheme. The central line density is p = (32 ± 7) atoms//xm (error 
given by standard deviation) in each well at an initial temperature of the unsplit 
cloud of T = (120 ± 30) nK. For the initial gas before the splitting this results in a 
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Figure 7. Comparison of distributions of the squared contrast obtained with 
the optical sHcing method and the direct imaging system for the same double- 
well configuration. The evolution time is 10ms and each distribution contains 300 
realizations, which minimizes statistical effects. In the upper row we plot exemplary 
distributions of the squared contrast, in the lower row the squared contrast is 
normalized to its mean. While the former shows slight differences for the results of 
both imaging systems, the latter shows very good agreement (see text for details). 

chemical potential of /i ~ {2/3)hur and an energy fc^T ~ hur due to temperature. 
These parameters are similar to the parameters of other experiments with Id gases in 
microtraps, e.g. [72, 73], and give an effectively Id Bose gas in the quasi-condensate 
regime. 

After the splitting, a slight net mean imbalance An = {Nl — Nji) / {Nl + N^) of 
An = (1.5 ± 1)% is present. Here, A^^^ and Nji are the atom numbers of the left and 
right well, respectively. Due to technical noise the width of this imbalance distribution is 
around 2 — 3 times larger than the quantum shot noise limit of a/ZV, where N = + Nji 
is the total atom number of the system. The imbalance has mainly two effects. The 
net difference in chemical potential resulting from the mean imbalance leads to a trivial 
common phase evolution of the two wells [37, 74, 75]. In addition to this, the width 
of the imbalance distribution leads to an additional broadening of the global phase 
distribution. It is a strength of our approach that this technical noise does not affect 
the dynamics of the contrast distributions, as its effect on the dynamics is negligible 
compared to the quantum noise associated with the splitting process for the probed 
integration lengths L. 

For the shortest evolution time te = 2 ms the C^/(C^) distributions are peaked 
on all lengths L, being almost all identical in form and having a low probability of 
obtaining a low contrast C{L). As time evolves, the distributions behave very differently 
on different length scales. For long integration lengths L, the FDFs evolve to become 



CONTENTS 



16 



(a) FULL CLOUD 60 |jm 
tE=2ms 



38 |jm 27 |jm 22 |jm 




(b) 



01230123012301230123 

Normalized Squared Contrast C^/<C^> 

20 



15 



T3 

o 

T3 





• 




• T=T 

best fit 

• T=30nK 

• T=120nK 

^red"'' 




® • 

• 














o 


• 

• 


• 


• 

• 

• 


• 

• 


o 
• 






• • — 





10 20 
Evolution Time (ms) 



30 



Figure 8. (a) Histograms for the distributfons of C^/(C^) for evolutfon times from 
2 to 27 ms after splitting, shown together with the respective best fitting equilibrium 
distributions, which are Toff = Tbcstfit = (14 ± 4), (17 ± 5), (14 ± 4) nK for te = 
12, 17, 27 ms, respectively. For comparison, the calculated equilibrium distributions 
for T = 120 nK (blue dashed line) and 30 nK (green dot-dash line) are added, (b) The 
evolution of the best fit's value for the reduced with hold time. The probabilities of 
observing values of the reduced chi^ which are larger than those calculated from the 
best fits of the measurements are 50%, 62%, 18% for the last three evolution times of 
te = 12, 17, 27 ms respectively. 
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exponential in form, i.e. there is a high probabihty of observing a low contrast C{L) and 
the initial correlations of the system appear to be lost. For short integration lengths, 
however, the FDFs remain peaked and there is a very low probability of observing a 
low contrast C{L). This directly demonstrates the persisting memory in the system 
of the initial state. Furthermore, the strong length dependence is a direct signature of 
the multimode nature of Id Bose gases [40, 76]. Finally, the shape of the distributions 
evolves very quickly for the first 12 ms and very little between te = 12 — 27 ms, showing 
that the system reaches a quasi-steady state. This imposes the question of whether this 
state corresponds to the thermal equilibrium state of the system. 

4-2. Emergence of an Effective Temperature 

We will now address the question of whether the observed quasi-steady state already 
corresponds to the thermal equilibrium state of the system. 

If the system has reached thermal equilibrium, the measured distributions plotted 
in figure 8a should be consistent with equilibrium theory [35, 56, 41]. The obtained non- 
equilibrium distributions of C^/ (C^) were therefore fitted with equilibrium distributions 
taking temperature as a free fit parameter. To this end, a analysis [77] was performed. 
Figure 8a shows the experimentally obtained distributions together with the respective 
best fitting equilibrium curves. The agreement is very good for longer evolution times, 
but less good for early evolution times. 

Figure 8b shows the corresponding reduced values of the fits. For short evolution 
times the reduced differ significantly from 1, showing that the experimentally 
obtained distributions are inconsistent with equilibrium theory. For evolution times 
longer than 12 ms, the observed value of the reduced approaches 1, which shows 
that the experimental data agrees well with equilibrium theory. This indicates that the 
thermal-like nature of the distributions is established dynamically during the evolution. 

It is thus possible to associate the observed steady state with a temperature T^g 
extracted from equilibrium theory. For the last three evolution times t^, = 12, 17, 27ms 
the fits give Tcs = (14±4), (17± 5), (14±4) nK, respectively. This is, however, a factor 
of eight lower than the initial temperature of the unsplit system. Furthermore, the 
corresponding effective thermal phase coherence length of the relative phase between 
the two gases is = p/mkBT^s = 13l|, lOtg, 13l|/im. This is much longer than 
the expected Xt = (1-5 ± 0.6) /^m deduced from the initial temperature of the unsplit 
cloud. For a detailed study of the emergence and properties of this length scale see [76] . 

As a further illustration of the striking difference between the observed FDFs and 
the ones expected for the corresponding equilibrium system, we performed interference 
experiments with two independently created Id gases, employing the same double-well 
potential as for the non-equilibrium situation. The measured equilibrium distributions 
are plotted in figure 9. Even for T = (60 ± 20) nK the measured C^/ (C^) distributions 
are exponential in form for all length scales probed, which is in stark contrast to the 
observed non-equilibrium steady state distributions but in perfect agreement with what 
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is expected from equilibrium theory. A second, colder dataset at T = (27 ± 7) nK 
visualizes the crossover from an exponential to a Gumbel-like shape of the equilibrium 
distributions [32, 56]. 
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Figure 9. Experimental (histograms) and tlieoretical (red solid lines) equilibrium 
distributions, (a) T = (60 ± 20) nK, density of (32 ± 6) atoms//ini. Each plot 
contains the statistics of ^ 170 measurements, (b) T — (27 ± 7) nK, a density of 
(27 ± 5)atoms//im, statistics of ~ 90 measurements. Both datasets show very good 
agreement with the theoretical description. 

It is, however, a justified question as to whether the very low observed T^s of the 
quasi-steady state could also be explained by a cooling effect involved in the splitting 
process or, stated differently, how does the initial temperature of the unsplit cloud relate 
to the final temperature of the gas after splitting and after a possible equilibration? In 
principle the rapid splitting procedure is a rather violent process which is very likely 
to transfer some energy to the system. This would suggest that the final equilibrium 
temperature should be even higher than the initial temperature of the unsplit system, 
in contrast to the steady-state that we observe. 

On the other hand, the sphtting process leads to a decompression of the gas, 
which, for a 3d ideal gas would lead to a temperature reduction. If we assume an 
adiabatic splitting a comparison of the trap frequencies leads to a lower bound on the 
decompressed temperature of (0.6 ± 0.1) x T [78]. This temperature reduction has 
been confirmed by preparing a thermal 3d cloud of atoms in the initial single-well trap 
and turning on the double- well potential using the same 17 ms ramp as in the non- 
equilibrium measurements. For this procedure, we observe a decrease in temperature of 
(0.59 ± 0.15) X T in agreement with decompression. 

To our knowledge, no quantitative descriptions of this process exist for degenerate 
Id gases. The exact evolution and distribution of the thermal energy during the 
splitting therefore remains a topic of ongoing research. The only quantities one can 
compare to are the initial temperature T = (120 ± 30) nK of the unsplit system and 
the lower bound of T ~ (70 ± 20) nK calculated under the assumption of an adiabatic 
decompression. Figure 8 therefore also shows the obtained values for the reduced 
if the experimentally obtained non-equilibrium distributions are fitted with equilibrium 
theory for T = 120 nK. The plot shows that these high temperatures are clearly rejected 
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by the test. In addition, in figure 8 the corresponding equihbrium distributions for 
T = 30 nK (green line) are plotted for comparison. This shows that even for T = 30 nK, 
which is much lower than the lower temperature bound obtained from decompression, 
the equilibrium distributions are still significantly different from what is observed in the 
experiment for the non-equilibrium system. 

These observations directly demonstrate that the quasi-steady state we observe 
has thermal-like properties but is not the true thermal equilibrium of the system. We 
thus associate this quasi-steady state with prethermalization, as introduced in [7] and 
suggested for split Id Bose gases in [24]. The system decays rapidly to a state whose 
contrast distributions are thermal-like in form, exhibiting a temperature almost an order 
of magnitude lower than the temperature of the initial unsplit system. As we will see in 
section 6, the relaxation of the system to this prethermalized state is very well described 
by a dephasing of the momentum modes of the system [23, 34, 76]. We will first briefly 
discuss the persistence of this prethermalized state as the system continues to evolve. 

5. Long-term Evolution of the System 

As shown in section 4, the non-equilibrium system decays rapidly to a prethermalized 
state. However, the question still remains as to any subsequent evolution of the system, 
particularly whether the system will eventually reach the true thermal equilibrium 
corresponding to two independent gases. 

In the measurements presented in the following, the initial temperature of the 
unsplit system was (78 ± 10) nK and the peak atomic hne density of the split system 
was (38 ±9) atoms/ /im in each of the two wells, i.e. approximately similar parameters to 
those presented in section 4 above. Furthermore, the final amplitude of the current for 
the dressing RF creating the double-well was increased to 23.8 mA leading to a higher 
potential barrier, larger separation of the double-well and therefore placing an even lower 
upper bound on any hypothetical residual tunnel-coupling. To ensure that the actual 
splitting process was also the same, the same linear ramp speed as in the measurements 
in section 4 was used, but the total ramp duration was extended to 18 ms. 

Figures 10 and 11 summarize the main results of the measurements. Figure 10a 
shows the evolution up to a time of almost 200 ms of the mean squared contrast (C^) for 
the full length of the system. The fast initial decay to the prethermalized state, followed 
by a very slow further evolution, is clearly visible. As an inset, example interference 
pictures obtained using the direct imaging system are shown. For very short evolution 
times the fringes are almost straight, directly visualizing the coherence of the splitting 
process. For larger evolution times the fringes become more and more wiggly, leading to 
the observed strong decrease in the integrated interference contrast. Again, for evolution 
times after the initial decay, the distributions of C^/ (C^) can be described by equilibrium 
theory, as shown for example in figure 10b for various evolution times te = 12 — 182 ms. 

Figure 11a displays the minimal reduced values for the inspected evolution times 
which shows that for evolution times after the initial rapid decay, the reduced is 
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Figure 10. (a) Time dependence of tlie mean squared contrast (C^) [8, 76] together 
witli example pictures obtained using the direct imaging system. (b) FDFs of 
C^/ (C^) extending to longer evolution times. Again, the experimentally obtained data 
(histograms) can be very well described by equilibrium theory where the distributions 
obtained at the best fitting temperature are plotted using red solid lines. The blue 
dashed lines show theoretical equilibrium distributions for the initial temperature 
T = 78 ± 10 nK of the unsplit cloud. 



similar to 1. The extracted effective temperatures are again very low compared to the 
initial temperature of the system, while higher temperatures clearly fail to describe the 
measured distributions, as demonstrated by the analysis for larger temperatures in 
figure 11a. 

Observing the long slow decay of (C^) seen in figure 11a one may be tempted 
to assume that this is the system slowly thermalizing. In order to analyze this slow 
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Figure 11. (a) Evolution of the reduced with hold time in the double- well 
potential, (b) Evolution of Tcff for the whole data set extracted using fits to theoretical 
equilibrium distributions (figure from [8]). The dashed linear fit indicates an increase 
of Teff over time of (0.14 ± 0.04)nK/ms. The colored area indicates the measured 
heating rate of the atom trap of 0.11 ± 0.06nK/ms, with the area being bounded by 
the error. 



subsequent evolution of the quasi-steady state, the effective temperatures for all times 
after the initial decay were extracted, as shown in figure lib. As can be seen, T^s 
rises slowly over time at a rate of (0.14 ± 0.04)nK/ms. This is, however, comparable 
to the measured heating rate of the atom trap of (0.11 ± 0.06)nK/ms which was 
characterized independently using equilibrium quasi-condensates. This indicates that 
either no thermalization is present, or, if it is present, that it is a very subtle process. 

6. Luttinger Liquid Description of the Dephasing of the System 

In the last sections it was shown that the system created through the coherent splitting 
of a single quasi-condensate rapidly relaxes to a quasi-steady state that has thermal 
characteristics but which does not correspond to the true thermal equilibrium state of 
the system. This observation can be very well understood by the theoretical model 
presented in [23, 24] where the appearance of a thermal- like state is explained through 
the dephasing of the multimode system. Before presenting a quantitative comparison 
of this model to our experiment we summarize the most important points of this theory 
in order to give an intuitive picture for an understanding of the prethermalized state. 
For a more detailed discussion the reader is referred to the original work [23, 24]. For 
an alternative description of the rapidly split quasi-condensate in terms of an effective- 
tunneling model, see [79] . 

6.1. Summary of the Theoretical Model 

The evolution of the interference pattern is directly determined by the evolution of 
the relative phase between the two uncoupled halves of the split gas. We describe 
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this evolution using a Tomonaga-Luttinger liquid approach. In this low-energy 
approximation, the relative degrees of freedom (/)(|/) and n{ii) perfectly decouple from the 
common degrees of freedom of the system 0com(?/) and ricom(2/), where the experimental 
observables are now described by operators. The resulting Hamiltonian for the relative 



where 



degrees of freedom is of the form H = ^ /_^i/2 '^^ ^i.^4^iy)Y + ^^"^iv) 
K = n^hp is the Luttinger parameter, c is the speed of sound, is the healing length 
and p is the atomic density. 

The evolution of (f){y) can be described in Fourier space by a set of decoupled 
harmonic oscillators of collective modes with momentum k. A collective mode 
with momentum k modulates the relative density and relative phase along the 
condensate in a sinusoidal fashion on a length scale of ~ 1/k. This is visualized in 
figure 12. 

In the experiment the splitting is performed fast in comparison to the timescale 
for the spread of perturbation which is set by the inverse chemical potential tgpiit < 
^h/c = Thus there is no time for the atoms to develop correlations along the 

longitudinal trap axis such that the relative density fluctuations (77,^)4^=0 = p/2 are 
completely random and not affected by atomic interactions. In particular, there is no 
correlation between modes with different momenta. Also, for each atom the decision of 
going to either half of the split system is random and uncorrelated with other atoms, 
leading to a binomial distribution of atom number fluctuations in each small segment of 
the Id system. The respective width of the relative phase distribution follows from the 
Heisenberg uncertainty relation. Initially, at te = 0, the relative phase fluctuations are 
governed completely by the atomic shot noise and are practically negligible at length 
scales > ^h- In other words, (0^)tg=o ~ and its fluctuations are strongly suppressed. 
On the contrary, initial density fluctuations exhibit for \k\ < ^ ^ large excess of noise 
compared to the zero-temperature state of two split quasi-condensates. In the latter, 
stationary case the relative density fluctuations are suppressed by atomic repulsion. As 
a result, almost all the energy is initially distributed in the density fluctuations, with 
each collective mode k containing equal energy but different populations that scale as 
l/\k\ as the energy per quanta in each mode is oc c|A;|. In particular, this leads to the 
weaker effects of modes with higher momenta on the evolution of the system. 

During the dynamics of the system the energy of each harmonic oscillator mode 
k oscillates between the fluctuations in density and fluctuations in phase, driven by 
interactions between the atoms. This results in a harmonic time dependence of the 
fluctuation amplitude with a period oc l/c|/c|. For short evolution times all phase 
fluctuation amplitudes grow in magnitude, which leads to a scrambling of the relative 
phase (f){y) along the longitudinal direction. This results in an initial rapid decrease in 
the interference contrast C{L). For longer times, the oscillations in different fc-modes 
dephase and the system reaches a quasi-steady state which is thermal in its appearance. 
In particular, correlation functions take a form that is algebraically equivalent to their 
respective counterparts in equilibrium. The effective temperature of the quasi-steady 
state corresponds to the energy which was equally introduced to all fc-modes by the 
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Figure 12. Illustration of the filter behavior of the choice of integration length during 
detection. If only one fc-mode was populated this would lead to the sinusoidal varying 
pattern displayed in the figures. Large fc-modes as shown on the left side affect the 
phase within the integration length in a uniform fashion almost not affecting the 
contrast. On the other hand high fc-modes as shown on the right produce a strong 
variation of the phase within L, strongly diminishing the contrast. 



splitting process. In analogy to the equipartition theorem [80] one can deduce that the 
quasi-steady state can be described by an effective temperature T^s given by [24] 

ksT,n = I = f . (4) 

where is the chemical potential, g is the Id interaction strength and fc^ is Boltzmann's 
constant. As shown in [24] the detailed form of the FDFs of C{L) and C^{L) can then 
be obtained by sampling the amplitude of the fluctuations of the phase with Gaussian 
statistics. The choice of a particular integration length L during the detection process 
represents a filter for the effects of different fc-modes (figure 12). For an integration 
length L, modes with wavelengths 2% /k < L produce strong fluctuations of within 
L, while modes with wavelengths 27r//c > L uniformly affect (j){y) in L and just change 
the overall phase of the interference pattern. This means that the contrast C of the 
integrated interference pattern is governed by the dynamics of modes with wavelength 
< L. 

The higher the population of a particular fc-mode, the greater the mean amplitude 
of its phase fluctuation. For long integration lengths L, there are many significantly 
populated modes with k > 2tt/L which affect the integrated interference contrast C. 
With increasing evolution time t^, the dynamics of these large amplitude modes leads 
to a strong decrease in the probability of observing a high contrast C and the FDF of 
the contrast takes the form of an exponential decay. This is the contrast decay regime. 
For short integration lengths L, the contrast C is influenced only by sparsely populated 
modes with high momenta so the FDF of the contrast preserves a non-zero peak for the 
whole evolution time. This is the phase diffusion regime. These two regimes have been 
experimentally studied in detail in [76]. The observation of these two regimes is a direct 
and intuitive visualization of the multimode nature of Id Bose gases [40, 76]. 

Note that the above discussion is only valid in the case of a perfect symmetric 
splitting. As discussed in detail in reference [24], a density imbalance between the left 
and right wells leads to a coupling between the relative and common degrees of freedom. 
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In this case, the initial temperature of the unspht cloud, which defines the initial state 
of the common degrees of freedom, influences the dynamics of the relative degrees of 
freedom and hence the dynamics of the interference pattern. However, the treatment 
in terms of collective excitations of the system remains valid along with the distinction 
between phase diffusion and contrast decay dynamics depending on the integration 
length L. The decay of the contrast then happens in two steps. First there is an initial 
rapid decay analogous to the symmetric splitting case. Secondly, the quasi-steady state 
which is reached, very slowly evolves in time due to the coupling between relative and 
common modes [24]. This is, however, not expected to lead to true thermalization as 
the system remains integrable and the modes will not fully equilibrate. 

6.2. Quantitative Comparison Between Experiment and Theory 

A quantitative comparison between experiment and the Tomonaga-Luttinger liquid 
model is presented in figure 13. We find good agreement of the theoretical description 
with the experimental data without any free fitting parameter. 

Considering that the theory neglects the effect of the longitudinal trapping potential 
and contains no fit parameter, only the experimentally measured input parameters, the 
agreement of experiment and theory is very good. For the first two evolution times 
the description is, however, less accurate for two reasons. First, the extraction of the 
contrast is very sensitive for the early evolution times, leading to a broadening of the 
experimental distributions by fitting uncertainties. Secondly, the disagreements also 
indicates that the initial preparation of the system through coherently splitting a single 
cloud is not fast enough to be instantaneous and perfectly uniform as assumed in theory. 

By preparing Id Bose gases with a wide range of different parameters, we can 
further systematically check the predictions of our theoretical model. Figure 14 shows 
the predicted scaling (Equation 4) of the effective temperature with density, as well as 
the independence of Tgfj from the initial temperature before splitting. 

These two observations provide further strong support for our theoretical model 
and the interpretation of the observations presented in this work as prethermalization. 
Especially the observation that the properties of the quasi-steady state are independent 
from the initial temperature of the unsplit system is another indicator that the observed 
state is clearly different from the thermal equilibrium state. The observed linear scaling 
of Tgff with density p further shows that the prethermalized state is solely defined by 
the quantum noise associated with the splitting process. The apparent small systematic 
offset of the experimentally obtained Tgfj and the theoretical description visible in figure 
14a can be attributed to imperfections in the experimental splitting process [81]. 

Returning to the discussion of the theory in section 6.1, one sees that the relative 
degrees of freedom and n do indeed decouple from the common degrees of freedom 
0com and ncom; at least on the timescales on which we observe the system. This is 
particularly interesting, since there is no a priori reason to assume that the different 
degrees of freedom should not couple and equilibrate in our system. 
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Figure 13. Comparison of distributions of the normalized squared contrast C^/(C^) 
obtained in section 4 to the theoretical description of the rapidly split quasi-condensate 
presented in the previous section 6.1. The experimental data is plotted using 
histograms and the theoretical simulations solid red lines. The light red shaded areas 
denote the errors resulting from the uncertainty of the experimentally measured theory 
input parameters. These input parameters were the experimentally measured values 
of density in a single well p — (32±4) atoms//im, trap frequency w_l — (1.4±0.1) kHz, 
imbalance An given in percent as (0.1 ± 0.7)%, temperature of the unsplit system of 
T ~ (78 ± 10) nK and the uncertainty of ±0.5 ms for the point in time of = 0. 
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Figure 14. a) Dependence of Toff on density p and b) Independence of Tog from 
the initial temperature T of the system before sphtting, corrected for the scahng of 
Toff with density. The (black) solid line corresponds to the theoretical prediction 
ksTcs = 9P/'^- The black (green) data point correspond to the dataset presented in 
section 4 (5) respectively. The figure was taken from [8]. 



7. Dephasing, not Thermalization 

In an earlier experiment by Hofferberth et al. Nature 449, 324 (2007) [38], the dynamics 
of a coherently split Id Bose gas was analyzed through the decay of the coherence 
factor \E'(t) = {Re J e'^^'^^'^^dz) . The coherence factor was first introduced in [34] and 
is closely related to the mean contrast. Note that while the observable C{L) is more 
easily extracted from experiments, it is more complicated to calculate theoretically. A 
characteristic sub-exponential decay of the coherence factor \l/(i:) oc exp[— (t/to)"] (with 
a ~ 2/3) was observed, as predicted by Burkov et al. [82] for the thermalization of such a 
system, where to is the characteristic timescale for the thermalization. This observation 
was interpreted as evidence for thermalization. 

The basic physical system investigated in [38] was very similar to that in the 
new experimental setup presented in this paper. However, at the time of [38], neither 
the experimental methods nor the theoretical understanding was sufficiently developed 
to make a full analysis and characterization of the many-body state to which the 
system decayed. As demonstrated in the present work, this characterization was now 
possible through the new, more complete theoretical description of the non-equilibrium 
processes [23, 24] and the ability to experimentally measure the FDFs which enabled us 
to characterize the state to which our system decays. We would like to revisit here this 
old work [38] in the light of our new results. 

In the work of Hofferberth et al. [38] the coherent splitting in was done in the 
vertical direction, along the direction of gravity. This allowed the direct observation of 
the phase fluctuations along the length of the system (as with the direct imaging system 
in this work, but without having to deal with the problem of the chip surface being in the 
line-of-sight of the imaging beam). However, in general, this simplifled imaging scenario 
comes at the price that a vertical splitting process is not very well controlled and non- 
symmetric due to the interplay between the pull of gravity and the 1/r dependence of 
the amplitude of the RF radiation emitted from the atom chip wires. As an illustration, 
a 1 /xm height difference corresponds to change in potential energy of more than 2 kHz, 
which is more than the chemical potentials in the Id gases. This gravitational sag can 
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be partially compensated by adjusting the polarization of the RF fields, by changing the 
distance of the trap to the chip surface and by taking advantage of the 1/r dependence 
of the trapping and RF fields used to generate the double- well potential [83]. However, 
this constrains the double-well splittings that can be implemented. Furthermore, the 
compensation is never perfect and variations of the potential depth during the spitting 
process remain that are on the order of the chemical potential of the Id gases. 

As a consequence, significant shot-to-shot variations in the atom number difference 
between the two wells were observed in the previous work [38]. These overall atom 
number fluctuations led to an additional phase evolution A0(t) = |A/i, where A/i 
is the difference in chemical potential between the two wells in each experimental 
run [37, 74, 75]. Consequently, in [38], these common fluctuations were rejected by 
evaluating the coherence factor through ^(t) = {Re J e^('t>{y,t)-Am)dy) = (| J e''^^y'^^dy\). 
Here A0(t) was determined for each individual interference pattern via A0(t) = 
arg (J d?/ e*'^^^'*'') , which for an inflnitely long system makes the expression for the 
coherence factor equivalent to the original formulation of Burkov et al. [82] . This way of 
analyzing the coherence factor, however, also rejected the contribution to the dynamics 
coming from the k = mode (i.e. the phase evolution which, in a semi-classical picture, 
can be interpreted as being induced by the quantum shot-noise of the difference between 
the total number of particles in each well). For a flnite size system this k = contribution 
to the phase evolution is suppressed only by 1/ where N is the total number of 
particles in both wells [34]. 

In addition, in [38], a strong longitudinal breathing mode was observed after the 
splitting which limited measurements to the flrst 11ms of the evolution. As discussed 
in section 3 we split horizontally, perpendicular to the direction of gravity in the new 
experimental setup presented here. In this way, the symmetry during the splitting 
process is easily conserved, which makes the initial non-equilibrium state much better 
deflned and the flnal evolution more repeatable. The atom chip in the new setup is 
designed in such a way that the longitudinal conflnement can be varied independently 
of the radial conflnement. This makes it possible to choose combinations that allow the 
splitting of a single trap into a double-well with a minimized excitation of longitudinal 
breathing modes. This made the probing of long evolution times possible. 

Furthermore, in the new experimental setup, in comparison to [38], we split much 
further away from the atom chip surface, facilitating a double-well separation that is 
much greater for a given barrier height than in the experimental setup detailed in [38] . 
Simulations suggest that this leads to a 10 times faster break down of the tunnel coupling 
between the wells during the splitting process because, due to the larger separation, the 
overlap of the wave-functions is much smaller once the classical connection between the 
wells breaks down. Hence, the requirement of fast symmetric splitting to reach the 
initial state which allows the observation of prethermalization (discussed in [23, 24]) is 
better met in this work. 

Any comparison to current theoretical models of the fast initial evolution requires 
a clear starting point of the evolution as the exact value of the characteristic exponent 
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a is very sensitive to the position of te = 0. By analyzing the new data and from novel 
insight from recent experiments on tunnel-coupled systems [33], it became clear that 
it is very difficult to determine the exact point in time when the two clouds separate 
during the splitting, the dynamics of the different momentum modes during the final 
moments of the splitting, and when an independent evolution starts. In the present 
evaluation we used similar criteria as in [38] to set the starting time of the evolution 
(simulated plasma frequency ^ 10ms and an observable phase drift). We find that the 
starting point of the evolution is very difficult to define to better than ± 0.5 ms, given 
our current experimental precision and theoretical understanding. This is not a problem 
for observing and analyzing the FDFs and the long-term behavior, but it prevents us 
from reliably evaluating a from our data. However, for completeness, taking the first 
9 ms of the evolution presented in section 5 and analyzing the coherence factor in the 
same manner as it was performed in [38], we extract a value for a that is compatible 
with the previous result. Yet, studying the system for different initial temperatures, we 
observe an independence of the characteristic relaxation time to from temperature, in 
stark contrast to the prediction of Burkov et al. [82]. For an ideal, sudden splitting, our 
model predicts an exponential decay with a = 1 [24, 34]. In the experiment the time for 
the splitting is finite. Thus, to settle the very interesting question of how the observed 
dynamics starts and to study its details will require extensive additional theoretical 
work on the dynamics of the splitting process itself, as well as advanced experimental 
techniques [84]. 

In view of the present experiment and the theoretical analysis, the interpretation of 
Hofferberth et al. Nature 449, 324 (2007) [38] should therefore be revised as showing the 
fast integrable dephasing of the Id systems, and not thermalization, as was suggested 
in [38] by the comparison to Burkov et al. [82]. This revised view is in line with the 
observation of the absence of thermalization in Id gases in optical lattices [6] . We point 
out that for the present experiment, even independent of our theoretical model, the 
observed clear difference between the FDFs of normalized squared contrast obtained 
for an equilibrium system formed by cooling in the double- well (figure 9) and the FDFs 
observed after dynamic splitting (figure 8) provides direct experimental evidence that 
no thermalization is observed in our experiment. 

8. Conclusion 

We have detailed the experimental observation of the non-equilibrium many-body 
phenomenon of prethermalization. This was accomplished by observing the evolution 
of a rapidly and coherently split one-dimensional Bose gas. Through the use of full 
quantum mechanical probability distribution functions (FDFs) of the contrast of matter 
wave interference patterns, the quasi-steady state to which the system decays was found 
not to be a thermal equilibrium state. Using theory developed for equilibrium systems, 
we were able to show that the FDFs of the interference contrast in the quasi-steady state 
are thermal-like in form, but display an effective temperature that is independent of the 
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initial temperature of the gas before the sphtting process. The early-time evolution of 
our system is well described by a Tomonaga-Luttinger liquid model that describes the 
dephasing of the system to the prethermalized state. This evolution can be visualized by 
means of FDFs and then contrasted to the true thermalization expected at long times. 
Moreover, these results clearly illustrate the power of this new method of probing the 
dynamics of many-body quantum systems through FDFs. Beyond the present work, 
FDFs are a very general and widely applicable method that can be extended to other 
quantum many-body systems [85]. 

The nature of the true thermalization that manifests itself as a slow evolution 
towards a steady state after the initial rapid decay to the prethermalized state still 
remains an open question. A possibility would be the integrability violation through 
virtual [86, 87, 88] or even real [86] excitations of the radial degrees of freedom. 

In fact, the question of thermalization in our system is closely connected to the 
question of how two quantum mechanically correlated but spatially separated objects 
can lose their memory of that initial correlations, i.e. how do classical properties emerge 
from a closed quantum system? To answer these questions of thermalization and the 
emergence of classicality is an ongoing experimental and theoretical endeavor. 
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